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Abstract. We present the direct measurements of electric-dipole moments for 
5^3/2 ~^ ^^5/2 transitions with 20 < n < 48 for Rubidium atoms. The measurements 
were performed in an uhracold sample via observation of the Autler-Townes splitting 
in a three-level ladder scheme, commonly used for 2-photon excitation of Rydberg 
states. To the best of our knowledge, this is the first systematic measurement of the 
electric dipole moments for transitions from low excited states of rubidium to Rydberg 
states. Due to its simplicity and versatility, this method can be easily extended to 
other transitions and other atomic species with little constraints. Good agreement of 
the experimental results with theory proves the reliability of the measurement method. 
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1. Introduction 

Rydberg atoms in the ultra-cold regime have attracted a lot of theoretical and 
experimental interest in recent years. The strong, long range nature of the interactions 
in these systems can be tailored with high precision, making them very attractive for 
applications in quantum information [U |2] and the study of strongly interacting many- 
body systems [2]. Crystal-like structures can be created where the atoms are confined in 
space and arranged in arbitrary geometries, and the interactions amongst them can be 
controUably switched on and off and tuned finely. The high level of control that can be 
achieved in these systems and the low level of decoherence they display has prompted 
several groups to study atom-atom entanglement and quantum gate operation using 
single atoms in microscopic dipole traps [2111]. 

The basis of qubit encoding and quantum logic gate operations is reliable 
preparation and manipulation of the electronic states of atoms and the tuning of 
their interactions. The recent demonstration of interaction-driven Rydberg blockade 
between two trapped ultra-cold atoms was based on the observation of an increase of 
the Rabi frequency in collective one-atom excitations [5]. Similarly, strong interactions 
in mesoscopic ensembles of ultracold Rydberg atoms can be exploited to create a novel 
state, where a number of atoms behave collectively to share a single excitation [6l [7]. 
The proof of the coherent nature of the excitation relies on the observation of Rabi 
oscillations between low excited and Rydberg levels. The Rabi oscillation of a single 
atom (or blockaded sample of atoms) between a low and a Rydberg state is governed 
by the dipole moment of the transition involved. In ref. p], Rabi oscillation of single 
atoms were numerically simulated, elucidating the coherent nature of the excitation and 
giving useful insight on the experimental limits. More generally, a comparison of the 
experimental measurements with theoretical calculations of dipole matrix elements is 
required, both for verification of the experiments and to constrain the limit of validity 
and the accuracy of present models. 

Since analytical expressions for radial matrix elements are known for hydrogen 
only [9], numerous theoretical models have been developed for the calculation of dipole 
matrix elements for alkali-metal atoms. Recent measurements of the effective lifetimes of 
Rydberg states [101 [H] can be used for indirect verification of the theory [I2] , but they 
do not reveal the spectroscopic features of each particular transition. Experimental 
measurement of the absolute transition probability between excited states of alkali- 
metal atoms is a challenging task [13]. The transitions between ground and first excited 
states have been studied in detail, both experimentally and theoretically. The oscillator 
strengths of the principal series of rubidium were measured [TU [15], as well as the 
oscillator strengths for transitions between low excited and Rydberg states of sodium 
and lithium, presented in ref. [Ml [17] . Most of the experimental methods have limited 
accuracy [T8j. 

More precise measurements of the dipole moments for transitions between excited 
states of alkali-metal atoms are therefore necessary for the verification of theoretical 
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Figure 1. Schematic representation of the three- level ladder scheme, ilc and ^Ip are 
the Rabi frequencies associated to the coupling laser and probe laser respectively. 



models. Particularly, measurements of the dipole moments for transitions between first 
excited and Rydberg states, which are important in many experiments with ultracold 
Rydberg atoms, have never been performed. 

In the present work we describe a relatively simple method, similar to one previously 
used in [19] for molecular transitions, which allowed us to measure for the first time the 
dipole moments for Rb 5P3/2 — ?■ transitions with 20 < n < 48. This method is 

based on the spectroscopic observation of the Autler-Townes splitting (AT) in a three- 
level ladder scheme [20], which is also commonly used in electro- magnetically induced 
transparency (EIT) experiments with cold Rydberg atoms [21], |22] . 

2. Method 

We observe EIT and AT splitting [23] in a sample of ultracold Rb Rydberg atoms by 
measurement of the absorption profile of the weak (probe) laser radiation scanned across 
the resonance 5Si/2 5P3/2 in the presence of the stronger (coupling) laser radiation 
tuned to transition 5P3/2 -> following the scheme suggested in [201 122] . 

The long-lived states |1) and |3) are coupled to a short lived state |2) by a weak 
probe laser and an intense coupling laser, as shown in figure [H and we assume that the 
states |1) and |3) are not coupled by electric dipole transition. With the atoms initially 
populating state |1) and the coupling laser tuned to the |2) to |3) resonance, EIT appears 
as a dip in the center of the absorption lineshape of the probe beam, which results from 
the destructive interference between excitation pathways, and is observed by sweeping 
the probe laser through the |1) — )■ |2) transition. With increasing coupling laser intensity 
the transparency dip becomes wider, until it eventually resembles two individual lines 
(AT splitting), resulting from the dynamic splitting of the middle energy level due to 
the strong coupling field. 

To obtain an expression for the absorption lineshape of the probe beam, we solve 
the optical Bloch equations for the system of three atomic levels and two coupling lasers, 
taking the probe laser as a weak perturbation so that the population of the the two states 
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1 2) and |3) are set to 0. This is ensured by keeping the probe laser intensity Ip <^ I sat 



where ctq is the cross-section for absorption in the absence of the couphng laser, F is 
the width of the |2) state, 5 the detuning of the probe beam, A the detuning of the 
coupling beam and 73 is the dephasing rate of the |3) state. In our experiment, state 
|1) corresponds to the 5S1/2F = 2,mp = 2 state, |2) to 5P3/2F = 3,mF = 3 and the 
state |3) to one of the Rydberg states nD^i2,F = A^rup = 4, where n lies between 20 
and 48. It is to be mentioned that other states (e.g. ^1^3/2) can be also selected, by 
tuning the frequency and changing the polarization of coupling and probe beam. For 
the 5^1/2^ = 2 ^ 5P3/2F = 3 transition, cxq = 2.90 x lO^^^ ^^-2^ 

By scanning the frequency of the probe beam around the |1) — t- |2) resonance, 
we obtain the spectrum o"p(5), from which we determine VLq- By investigating the 
dependence of Vtc on the power of the coupling beam for a variety of Rydberg states, 
we experimentally determine their transition dipole moments. We find that this is a 
highly reliable and precise method of determining dipole moments. As it will be shown 
in the next section, the method allows for a strongly focused coupling beam to be used, 
such that Vtc ~^ {x, y), provided an accurate determination of the laser beam profile 
and power is performed. 

3. Experimental details 

We prepare an ensemble of cold atoms at a density of 6 x 10^ cm~'^ and a temperature 
of 10 /iK, by collecting 3 x 10® atoms from the background vapour in a magneto optical 
trap and allowing further cooling in optical molasses. The sample is prepared in the 
|5S'i/2,-F = 2,mp = 2) magnetic sub-level by optical pumping, with an efficiency of 



The atoms are subsequently illuminated by an intense coupling laser and a weak, 
counterpropagating probe laser, as shown in figure [2^), which are a^/cr^ polarized. The 
coupling laser is provided by a frequency doubled diode laser (TOPTICA SHG) which 
is frequency locked to the center of the 5P3/2 — nDj transition at around 480 nm 
using the scheme based on KIT described in[2l]. The beam is focused onto the atom 
cloud to a small waist in order to maximize its signature on the atomic sample. The 
weak 'probe' beam is provided by a commercial diode laser frequency stabilized, which 
is imaged onto the cold atom cloud by a singlet lens (LI) to a waist of 364 ± 10/^m. 
Its frequency can be swept around the 55*1/2 — ^ 5P3/2 transition at 780.24 nm. After 
passing through the atoms, a system of two lenses formed an image of the probe beam 
on either a CCD camera or a photodiode, as shown in figure [2^). We are able to address 
nD'ij2 and nD^j2 states by tuning the laser appropriately, since the < IMHz linewidth 



for the transition to |5P3/2)- We find that the cross-section for absorption of the probe 
laser, ap is then: 




(1) 



80%. 




Figure 2. a) Scheme of the experimental apparatus showing the MOT vacuum 
chamber and the probe and coupling beams arrangement, b) Diagram of the tuning 
of the probe (red) and coupling laser (blue), showing the |5S'i/2, |5-P3/2, 1*^-05/2 states 
and the hyperfine structure (including magnetic sublevels) 



of our lasers is much smaller than the fine-structure splitting for the range of n in this 
study. However, we cannot resolve the hyperfine structure of these states. Thus, in the 
general case, we do not have a true three level system. By choosing both laser beams to 
have 0"^ polarization, we drive transitions such that Amp=l, according to the scheme 
in figure [2b). As the atoms initially populate the 155*1/2, -F = 2,mF = 2), the choice of 
polarization of the probe selects the |5P3/2,-F = 3,mp = 3) state. For the nD^i2 state 
this scheme selects only the F=4, m.p=4: sublevel. 

The atoms are illuminated for 1 ms by switching on both the 480 nm coupling and 
780 nm probe laser. During this time the probe laser frequency is swept across the 
551/2-^ = 2 to 5P3/2F = 3 transition, whilst a 200 mG magnetic field is kept on in order 
to preserve the projection of the atoms' magnetic moment. Finally, we record the total 
laser power for each realization of the experiment on a calibrated photodiode, accurate 
to better than 5%, also accounting for the response of the photodiode due to wavelength 
variation. We correct for the variations of the probe laser intensity while the frequency 
is swept. In this way we determine the probe absorption cross section as a function of 
its frequency using Beer's law. 

It must be pointed out that the Rabi frequency flc that appears in equation 1, is 
in general a function of position Qq ~^ u)- In many experiments, a homogeneous 

Rabi frequency across the sample is necessary and this requires proper shaping of the 
coupling laser beam. The analysis method that we employ allows us to avoid the 
reshaping of the beam, thus enabling us to maximize the intensities involved by focusing 
the laser beams. 

Our method, therefore, does not allow us to resolve the distribution flci^, u) of the 
coupling Rabi frequency due to the Gaussian intensity profile. Instead it provides an 



Rydberg transitions electric dipole moments 



6 




integrated signal and an 'average' cross section <Jp{6) across the spatial profile. Taking 
this into account in the analysis of the spectra is non-trivial as an analytic solution to 
equation 1 cannot be found. Thus the signal is numerically evaluated by considering 
the contributions from areas of equal intensity and integrating over the beam profile, 
see below. By taking into account the spatial profile E{x, y) and therefore fic(a^5 u)^ we 
extract from each scan a value for il™"^ which correspond to the peak value of E{x, y), 
i.e. the amplitude of the Gaussian distribution. Figure [3] shows how the modelling of 
the lineshape well adapts to the measured spectra for n=22 and n=44. 

For each chosen Rydberg state we take scans at different coupling powers, 

while the probe laser is kept at very low intensity (typically O.OlIsat) to minimize the 
probability of populating the 5P3/2 state, as mentioned earlier. This allows us to build 
a linear trend for each n of versus y/P. From the measurements we can therefore 
extract a value for the transition dipole moment for each n, as ^Ic oc /i \/P 

4. Analysis 

We have analyzed the probe absorption lineshapes to determine the Rabi frequency 
Qmax q£ ^YiQ upper transition. For an accurate analysis, it is fundamental for us to 
account for the inhomogeneity of the intensity profile of the coupling laser beam. A 
useful feature of our setup is that we may take absorption images of the atoms in 
the presence of the coupling beam, which allows for an accurate determination of the 
blue laser waist at the position of the atoms. The atomic cloud casts a shadow on 
the probe beam that is imaged on the camera. The presence of the coupling beam 
affects the optical depth of the atomic cloud, so that the cloud becomes transparent to 
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Figure 4. Snap shot of the rubidium atoms absorption imaged using the probe and 
the counterpropagating couphng laser. The hole in the center of the cloud is due to 
the reduction in cross section for the probe laser where sufficient coupling intensity 
exists to produce EIT. An accurate measurement of the probe laser waist is obtained 
by fitting the expected 2D profile to the observed image. 

the probe beam. By imaging the center of the atomic cloud, the optical depth of the 
atoms in the region illuminated by the coupling beam should reflect its intensity profile. 
Figure H] shows a typical variation in the intensity of the light transmitted through the 
sample. This allowed us to accurately determine the elliptical coupling laser waists as 
Wmaj = (240±10)/im Wmin = (172±10)/im, along the major and minor axes of the beam 
cross section. These measurements are mainly limited by the resolution and signal to 
noise ratio of our imaging system. To make sure that the measurement were not biased, 
we repeated the measurements for different coupling laser powers (80 mW and 50 mW) 
and tuning to different n states (n=26, 40, 44). We derived consistent waists in all 
cases, thus confirming the effectiveness of the method. Using this determination of the 
laser beam profile, and the measured total power of the beam, we deduced the electric 
field amplitude E(x, y) of the coupling laser as a 2-D Gaussian distribution of amphtude 
Emax and waists Wmaj and Wmin, which is then taken into account via numerical methods 
in the analysis of the spectra. 

In order to obtain good AT spectra, we were obliged to work with very low probe 
intensity, and consequently we compromised the bandwidth of the photodiode amplifier 
to 35 kHz in exchange for a higher signal to noise ratio. As a result, our spectra are 
instrumentally broadened and we are unable to observe the narrow EIT features as seen 
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in We quantified the broadening by analysing the spectrum obtained when the 

couphng laser is absent i.e. by observing the two level |5S'i/2-F = 2) — )■ |5P3/2,-F = 3) 
transition. We find that the linewidth F is 9 MHz, larger than the value of 6.065 
MHz obtained in high precision measurements. To explain the observed broadening, we 
modelled the experimental signal by convolving the ideal spectrum with the response 
of a passive low pass filter with a corner frequency of 35 KHz. We found that the 
photodiode bandwidth, combined with the measured laser linewidth of 450 kHz, was 
sufficient to explain the signal broadening. This instrumental limitation determines the 
width of the EIT feature that we can observe, and limits our measurements of 73 to 2.5 
MHz, much larger than the values expected from the lifetime of any the Rydberg states 
used. In order to check that the instrumental broadening did not alter the value of Qc 
obtained from fitting the lineshape, we studied the effect of the filter on the AT signal 
(taking into account the inhomogeneity of the coupling laser). We found that the our 
fitting procedure yielded the correct to high accuracy (better than 1%). 

We analysed the spectra as follows. In general, the photodiode signal is the given 
by Beers law as the integral over the atomic cloud: 



where the cross-section ap{x,y) (equation 1) is a function of x and y due to the 
inhomogeneity of the coupling laser and the atomic density n{x,y, z). We approximate 
this integral as a sum using as input the measured waist sizes of the coupling beam and 
the measured cloud dimensions. The spectra are fitted to this sum and we obtain the 
parameters and 73 and A, having already determined F as above. Figure 4 shows 
examples of the data and the fitted lineshapes. 

We repeated the measurements for each n state at different coupling laser powers, 
ranging from 5 mW to 80 mW. The results obtained for n = 22 and n = 44 are plotted 
in figure [5] together with their uncertainty, which is given from the statistics of the fitting 
procedure, and show a good linear trend. The Rabi frequency is in fact directly linked 
to the square root of the laser power, according to f2„ = /i-f = /i„ — . ^ \fP, where 

E is the amplitude of the electric field in a Gaussian beam of power P and waist w. In 
our case, because of the ellipticity of the beam, w"^ = Wmaj x "U^min- 

The dipole moments for each transition |5P3/2,-F = 3,mF = 3) — )■ \nDr,/2,F = 
4,mF = 4) (with n=22 to 44) can be therefore obtained directly from the gradients of 
the linear trends in graph [5l provided an accurate measurement of w and P is made. 
From the linear fits we obtained the values for the gradients with an uncertainty of less 
than 2% in most cases, with a maximum of 4% for n=44. However, our uncertainties on 
the measurements of the waist and of the laser power were also to be taken into account 
when determining the values of the dipole moment /i„ for the transitions involved, and 
the total uncertainty on the measured dipole moments results to be of less than 10% 
for all n. 

In figure [6] our results are plotted, after being rescaled, as reduced matrix elements 
to be representative of the I5P3/2) — )■ \nD5/2) transition, by taking into account the 




(2) 



Rydberg transitions electric dipole moments 



9 



N 



O 



25; 

20 
15 
10 



^ n=22 
I n=44 



S 5 







0.00 0.05 0.10 0.15 0.20 0.25 
^Coupling Power (W^^^) 



Figure 5. The dependence of the Rabi frequency ft on the square root of the couphng 
power VP is used to check; the inhomogeneity of the coupling laser has been correctly 
accounted for. Error bars represent the uncertainty in the value of f2 obtained from 
the fit. 
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Figure 6. Our results, summarized as a graph of measured reduced matrix elements 
fi versus principle quantum number n. Good agreement with 3 of the 4 models is 
obtained, however the method of Marinescu yields slightly higher values for /z 
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three-J and six-J symbols to obtain an angular coefficient of 




. The plot shows the 



reduced matrix elements and their uncertainties versus n-number. 

To demonstrate the validity of our measurement method, we compared our results 
to theoretical models. Many of the methods currently used for calculations of transition 
dipole moments include a relatively simple numeric integration of Schrodinger equation 
in the Coulomb approximation. Quantum defects are taken as input parameters or 
model potentials are introduced to account for the interaction of the Rydberg 
electron with the atomic core. One of the most popular methods uses the Coulomb 
approximation with the quantum defects of the Rydberg states taken as input 
parameters. The Shrodinger equation can be solved numerically [26] (Numeric Coulomb 
Approximation, NCA). The radial matrix elements can also be obtained analytically 
from the quasiclassical approximation [271 EH] (Kaulakys or Dyachkov-Pankratov DP 
model) or by extrapolation of the exact radial integral for hydrogen to non-integer 
quantum numbers [29] (Modified Coulomb Approximation, MCA). A second group of 
methods is based on the use of the realistic model potentials (Method of model potential, 
MMP), which takes into account the penetration of the Rydberg electron into the 
atomic core, the polarization and perturbation of the core. More complicated relativistic 
methods, based on Dirac equation with corrections to core polarization were used for 
calculation of oscillator strengths for alkali-metal principal series and transitions between 
excited states with n < 11 [301 EI]- These are the most accurate and the anomalies 
in the ratios of experimentally measured oscillator strengths in potassium, rubidium 
and cesium were correctly described [30]. Comparison of the calculated reduced dipole 
matrix elements [32] with the experiment is shown in Fig. 4. We calculated the dipole 
moments using DP , NCA , MCA and MMP methods, described above, using the quoted 
values for quantum defects for Rubidium from [33], and the energy of the I5P3/2) state 
from [34] and the model potential of Marinescu [25] with corrections for fine structure 
and core perturbation. For MMP calculations the RADIAL program, discussed in [35] , 
was used. We used a 'chi-squared' statistical parameter to find which model describes 
the experiment better. The best agreement between experiment and theory is observed 
for MCA model, but agreement with DP, NCA is also evident. Our results highlight 
a discrepancy with MMP calculations, assuming that the inherent limit in accuracy of 
this model is better than 5%, as quoted in previous works benchmarking the model. 
New measurements of the dipole matrix elements for other series in alkali-metal atoms 
will be of great importance for further verification of the theory. 

5. Outlook 

We performed direct measurements of the dipole moments for transitions from the 5P3/2 
state to high Rydberg state uDj for ^^Rb using EIT/AT in a three- level ladder excitation 
scheme. These are, to our knowledge, the ffist direct measurements of the alkali- atom 
Rydberg state dipole moments using this method. Good agreement 
between experiment and theory is observed, and the measurements could be extended 
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to different atomic levels and species, so that this method could prove invaluable for 
testing current theoretical models with limited accuracy. 
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